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Abstract 

"en; 

In the intersecting braneworld models, higher curvature corrections to the Einstein action are 
necessary to provide a non-trivial geometry (brane tension) at the brane junctions. By introducing 
such terms in a Gauss-Bonnet form, we give an effective description of localized gravity on the singular 
delta-function branes. There exists a non- vanishing brane tension at the four-dimensional brane 
intersection of two 4-branes. Importantly, we give explicit expressions of the graviton propagator and 
show that the Randall-Sundrum single-brane model with a Gauss-Bonnet term in the bulk correctly 
gives a massless graviton on the brane as for the RS model. We explore some crucial features of 
completely localized gravity in the solitonic braneworld solutions obtained with a choice (£ = 1) of 
solutions. The no-go theorem known for Einstein's theory may not apply to the £ = 1 solution. 
As complementary discussions, we provide an effective description of the power-law corrections to 
Newtonian gravity on the branes or at the common intersection thereof. 
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Qh! 1 Introduction 

In recent years, it has been understood that the brane-world models of non-compact extra dimen- 
i sions [1] reproduce the usual Einstein gravity. This idea was quickly generalized to a discrete patch 
of AdS±+N space accommodating N intersecting (AT + 2) branes [2]. Later, the RS singular 3-brane 
model with perturbations of bulk geometry was fully investigated in [3, 4, 5], and interesting discussion 
of the 'crystal' multi-braneworld solutions appeared in [6]. A class of solitonic brane solutions arising 
from a special fine tuning between the bulk cosmological term and the Gauss-Bonnet coupling was 
given in [7]. The behavior of gravity localized on branes or string-like defects in N > 1 dimensions was 
studied in [8]. These constructions have been useful in providing some realizations of trapped four- 
dimensional gravity to a brane embedded in an AdS space. It could be possible to localize ordinary 
matter and gauge fields to this (3 + 1) dimensional manifold if the RS brane is a Z?3-brane-like object 
of string/M theory. Since a D3-brane in the string theory background is necessarily a solitonic object, 
and only the higher curvature corrections can provide a non-trivial background at the brane junction, 
we find it interesting to study intersecting brane-world configurations by including a Gauss-Bonnet 
term. 



In the intersecting brane configurations of [2, 6, 9, 5], a 3-brane given by the four-dimensional 
intersection of (2 + iV)-branes (N > 2) with larger co-dimensions has been viewed as our universe. 
In this picture, though each extra dimension has an infinite extension, the effective size of N non- 
compact extra dimensions is supposed to be finite, which could be as large as 1/xm (for N > 2) or 
much larger than their natural scale £ >> Mp\ (~ lCT 33 cm) [10, 2]. While taking a Gauss-Bonnet 
coupling into account, the effective four-dimensional Planck mass reads Af?// = ^* +N Vn, where the 
effective volume of N extra dimensions is 



Here £ is a coupling constant arising from a GB term, and the constant k has dimension of inverse AdS 
length. One of the motivations here to include a GB coupling is that the low energy version of certain 
string theories (such as type I heterotic strings) admits the quadratic curvature corrections as a GB 
invariant in their effective actions [11]. It is remarkable that the inclusion of a GB term can consistently 
explain the basic features of the Randall- Sundrum type brane- world models in D = 5 [12, 13, 14]. 
Some other aspects of Randall-Sundrum brane models supplemented with higher curvature terms were 
studied in recent papers, including [15, 16, 17, 18, 19, 20, 21]. Here we construct brane-world solutions 
for N > 1 with a Gauss-Bonnet term. Indeed, the inclusion of a GB term in the bulk might shed 
more light towards the understanding of the RS models in the AdS^+N spaces. One of the interesting 
observations is that a GB term does allow a non-trivial four-dimensional brane tension when there 
are two extra transverse directions. Also a Gauss-Bonnet term could be useful to circumvent some of 
the arguments given for the no-go theorem [22, 23] applicable to the braneworld models based on the 
Einstein action. 

The rest of paper is organized as follows. In section 2, we give a general set-up of the effective 
braneworld action, and present a framework for intersecting braneworld configurations with a GB 
term. Section 3 explains the basic features of the linear perturbation equations in (4 + N) dimensions. 
In section 4, we derive a relation between the effective four-dimensional Planck mass and (4 + N) 
dimensional fundamental mass scale. In section 5, we give the expressions for graviton propagators 
in AdSi+N space, and discuss some salient features of gravity localized on the branes. In section 6, 
we elucidate upon the braneworld solutions in five dimensions, which exhibit that £ = 1 is indeed a 
physical choice. We then give in section 7 a complementary description of Newton potential corrections 
for solitonic branes. In section 8, we briefly comment upon the scalar and vector modes of the metric 
fluctuations. Section 9 contains conclusion. 

2 Intersecting branes setup 

In the low energy effective action of the Horova-Witten type I heterotic string theory [24], there arise 
additional interactions of the gauge fields and higher powers of the curvatures. In particular, there is 
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(1) 



R 2 interaction in the Gauss-Bonnet combination [25]. In many cases a contribution of higher curvature 
terms is neglected as it is higher order in derivatives, but in the warped braneworld scenario in D > 5, 
we see the relevance of these terms if they are expressed in a Gauss-Bonnet form. We will not stick here 
to the physical origin of a Gauss-Bonnet term. Rather, we will try to explore the physical importance 
of this term in warped backgrounds. We shall start with the following braneworld effective action in 
D{= (3 + AT) + l)-dimensional space-time (£>), where dB represents the (3 + AQ-dimensional boundary, 

S = j^ +N x - 2 A + a (R 2 - AR PQ R p Q + RmnpqR MNPQ )^ 

N pi'th brane 

+ W d 3 + N x^g^ (-A,) 

~1 JdB 



i=l 



^ rij'th brane 

+ E/ ' d 2+N x^g^ {-Aij 

i+3 



+ / d 4 xV-^ 1=0 '-' 2 " =0 (-A) . (2) 



Here P,Q, - ■ ■ = 0, 1, ■ • • , 3 + N, x M (/i = 0, 3) parametrize brane coordinates with Lorentzian signa- 
tures and z l (i = 1, N) parametrize extra non-compact dimensions. The action (2) characterizes an 
array of A" orthogonal (AT+2)-spatial dimensional branes in (4+N) dimensions. A is the D-dimensional 
bulk cosmological term, A; are the brane tensions of the A" intersecting (A r +2)-branes, Ajj are (N+l)- 
brane tensions, and A is the four-dimensional brane tension at the common intersection of higher 
dimensional branes. The coupling a has the mass dimension of , and k = 167tG4 + at = M 2+N . 

The graviton equations of motion derived by varying the action (2) with respect to g MN take the 
following form 

N . 

\/~9i+N (G M N + kH M N + nAg M N) = ~ » X K \] -g^+lP S ( z i) ^M$N9pq = ^ 



2 

1=1 



N 



K A 



y/-g 4 (w~<*»=<» 5{ Zl )5{z 2 ) ■ ■ ■ 5{z N )5» M 5» N g z ^-^=° , (3) 
where Hmn, an analogue of the Einstein tensor stemmed from the GB term, reads 



Hmn = -^9mn{R 2 -AR PQ R PQ + Rpqr S R PQRS ) 

+2a [RRmn - 2R M pnqR PQ + RmpqrRn PQR - 2R M P Rnp] ■ (4) 

In equation (3), the indices (p, q) take (3 + N) possible values, while (r, s) take only (2 + N) possible 
values, and fi, v = 0, 1, 2, 3. The (2 + AQ-brane tensions Aj multiply only S(zi) (i.e. the terms involving 
only one delta function), Aij multiply two delta functions, and the brane tension at the common 
intersection (A) multiplies a product of delta functions 6(z\), • • • , 5(zn). In D = 7, for example, A 



multiplies 5{zi)8{z2)5{z^) . One observes that when N = 2, is replaced by A, so that A can have 
a non-zero value at the intersection of two 4-branes as we exhibit below. 
We write a (3 + N) + 1-dimensional metric as 

ds 2 = e~ 2A{z) (g^{x x )dx»dx u + g ij {z)dz i dz j ^j . (5) 

We assume that the four-dimensional spacetime is Minkowski. Then we may express the metric (5) 
in the following form, which represents the Poincare half parametrization of AdS^+N space [2], 

t 2 

ds\ +N = ^2 (ji^ydx^ dx v + dz 2 + dfi^-i) . (6) 

The length scale £ is fixed by the (4 + AT) -dimensional bulk cosmological term A. From the exact 
non-linear analysis in (4 + AQ-dimensions we can find the most general solution to modified Einstein 
field equations, admitting a normalisable 4d graviton. This solution is A(z) = log(fej|zj| + 1), where 
Zi = Zi, Z2, ■ ■ ■ count extra spaces and ki = k±, k 2 , • • • are the inverse of AdS curvature radii. The 
choice k± = &2 = • ■ ■ = k further keeps the metric (6) manifestly symmetric under permutations of all 
extra N dimensions, and the length scale L (= (Vn) 1 ^) can be interpreted as the compactification 
size of N extra dimensions [2]. The metric solution then becomes 

1 N 

( k E i= i \*i\ + 1 ) 2 v j= i ' 

Here k = (yjnl)~ l . Since the coordinates parametrize the extra dimensions, the full bulk space 
[z? 0) comprises 2^ identical patches of the AdS^+N space, which are glued together along the 
branes (z 3 = 0). Here we have taken AdS space as the simplest choice of the bulk spacetime. And 
branes are characterized by the delta- function- like singularities, which arise from the second derivative 
of the warp factor A{z). However, there might exist other possibilities of an AdS bulk, such as the AdS 
Schwarzschild black-hole spacetime considered in [26]. It is also plausible that within the AdS/CFT 
correspondence gravity can be trapped on the brane in a manner similar to that in the original RS 
scenario [1], when the brane is embedded into the AdS Schwarzschild bulk [27]. It would be interesting 
to study the gravitational perturbations in such a black-hole background. 

We have assumed that the warp factor A(z) is a function of all transverse coordinates, that is, z % 



a{z) = \ og (jr>N +1 ) 



(8) 



Then a straightforward calculation yields 

N 

1/2 



A' 2 = e- 2A ^k 2 J2(dzM 2 = e- 2A(z) k 2 N, (9) 

i 

N 

A" = - e - 2A ^k 2 N + e- A( ^2kY,Hzi), (10) 



where primes represent differentiation with respect to z. The metric solution, in terms of the warp 
factor A(z), therefore satisfies a relation A" + A' 2 = e~ A ^ 2k5(z) > 0. Here, we simply note that 
this condition is equivalent to the weak energy condition one derives from G\ — G z z = k (t£ — T'f) = 
(N + 2) (A" + A' 2 ^j > 0. 

In order to exhibit localized behavior of gravity to the brane intersections, one may look at the 
linear perturbations about the background given by (7). For this purpose, we perform a conformal 
transformation gMN = e~ 2A ^gMN, where A(z) = log(fc^j \z{\ + 1). Then, we adopt a background 
subtraction technique as in [9]. Being specific, we use the linear equation 5Gmn + k5Hmn = 0, 
where SGmn = SGmn — Guph^ F ', 5Hmn = 5Hmn — HMph^ p . Note that Gmn = kTmn gives 
the background Einstein equations. What we are equating here with Gmn are only the vacuum 
energy contributions to the energy-momentum tensor. In our approach, the Randall-Sundrum-type 
fine tunings will only be implicit. If we wish, we can see them by considering the difference between 
(<5G>i, + k5H^ v ) and (5G^ V + k5H^ v ). This gives after some simplifications, 



1 r 
2 



Rh 



2Ru Kx 



an 
~2~ 



{R 2 



4Rmn 



R 



mnpq) 



4 (RR„ a - 2Rp pxq Rpq - 2i^ 5 i? AS + R^pqr R XPQR ) h uX 



(11) 



One could fine-tune the contributions coming from these terms to the vacuum energy contributions 
on the branes, that is, the source terms on the right-hand side of equation (3). For this purpose, we 
parameterize h^ u = e^ N+2 ^ A ^ z ^ 2 and obtain 

N 

N(N + 2)(N + 3) k 2 e~ 2A (1 - £/2) - 8(N + 2)k e~ A {l - f /3) ^ 5{ Zi ) 

h^,i>{x, z) = 



i=l 



-16(JV + l)(N + 2)aKk 2 J25{zi)5( Zj ) 

N 



-2k 



-2A(z) A 



+ e 



-A(z) 



i=l 



h/j,u{x, z) , 



(12) 



where £ = 2N 2 (N + 1) OL K M~ 2 k 2 and we have replaced a by M^a', so that a' is a dimensionless 
Gauss-Bonnet coupling. Here, it might be relevant to note that in the full non-linear analysis the 
fine-tuned relations emerge as boundary conditions to be satisfied on the brane(s), while the bulk part 
amounts to express A in terms of k 2 ; see [17] for the N = 2 case. For arbitrary TV, we obtain 

N(N + 2)(N + 3)fc 2 (1-^/2) 



A 



A, 



2k 

4(JV + 2)fc(l-£/3) 



A 



8(N + 1)(N + 2)a' Ml k 2 



(13) 

(14) 
(15) 



For N = 2, one may replace the brane tension Ai^ by A, because there is only one four-dimensional 
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brane intersection for two 4-branes that mutually intersect. A clear message is that the brane tension 
Ai ; 2 is non- vanishing for a' > 0, but it vanishes for the Einstein gravity (a = 0). 
We may perturb the background metric as 

ds l+N = n ^ N ] | 7T 2 { + V) dxV + E ( d ^) 2 | • ( 16 ) 

(*E<=i N + 1 ) 2 \ j=i / 

In this gauge, however, there may be some additional gravitational degrees of freedom coming from 
off-diagonal components (h^ m ) and diagonal components (h mm ) of the perturbed equations for Iimn, 
where m could run from z\ to zn- We will comment upon these modes in the section 8. One may 
also use the approximation where the tensor modes either decouple from the perturbed equations for 
h^m and h mm , or only the non- vanishing components of the fluctuations are h^ v [2, 5]. In fact, in 
the ordinary two-derivative gravity supplemented by a GB term there are no extra degrees of freedom 
than that of Einstein's theory. 

We find reasonable to study first the linear tensor fluctuations, h^ v in the gauge ^ = 0, d x h\ ll = 
0. The linearized equations for h^ v then take the following form 

N 

h„ 



i=i 



K 



+2aN 



N 



N(N + l)k 2 (D 4 + D z ) - Ak e A{z) S( Zi )a 4 - 4k e A{z) S( Zi )d 2 z 

i=l i^j 
N N 



+ (iV+ l)k 2 ^Y,Hzi) ~ N{N + 2)ke- A{ - z ^Ysgn{z j )d z 

i=l j=l 



V = ( 17 ) 



where U z = d zi + d Z2 + ■ ■ ■ d 2 N . When a = 0, our results will reproduce the results appeared in [2]. 
But here we find more interesting new features which are unavailable if the gravity action does not 
contain higher curvature corrections. 

3 Localized gravity in (4 + iV)-dimensions 

One can remove from equation (17) the single (linear in) derivative term from first square bracket 
by re-scaling the metric: h^v = e^ N+2 ^ A ^ Zi ^ 2 h^. We should note, however, that single derivatives of 
h^ v , and terms like 5(z) sgn{z)d z still survive from the second square bracket in (17). Fortunately, 
suitable combinations of these terms satisfy the necessary jump condition(s) on the brane(s). In fact, 
such terms do not appear if one has started only with the Einstein action, i.e. a = 0. We can make 
a change of variables \i^ v {x^z) = h^^x^ip^z) = t^ v ^ x ip(z}, where e^, is the constant polarization 
tensor of the graviton wave function, and p 2 = —q 2 = m 2 . Then, to the leading order, 

S~ j d N z \i/>(z)\ 2 J d 4 xd% l/ (x)d x h^{x) + --- . (18) 



This actually implies that (i) ip(z) is the conventional quantum-mechanical wavefunction, and (ii) if we 
had not chosen the coordinates in a conformally flat form, the kinetic terms in the x and z directions 
would have had different conformal factors, but the redefinition of the metric = e ( N + 2 ) A ( z i)/ 2 
puts the action in the canonical form by absorbing the conformal factors. This is similar to what 
happens in a model without the Gauss-Bonnet interaction term [5]. 
The metric fluctuations in terms of h^ v take the following form 



1 

K 



h 



+Nae 2A ^ 



2N(N - l)k 2 e~ 2A ^ (D 4 + - 8ke~ A ^Y, 6 ( z i) a 4 



+8 k 2 e~ 2A ^ ((N + 1) ^ sgniziWz^ - £ S(zj)sgn(zi)d g 



-8ke' A ^ £ 5{ Zi )dl. - 4(iV + 2) JbV^M k 2 e -2A( z ) 



+2^^)^)" (^ 2 + 2) e A ^Y^S(zi) 

i+3 i / 



h flu = 0. 



(19) 



Here we have defined n^h^ = m 2 h^ u . In terms of tp(z), we may express the above equation in the 
form 



N 



— m 



N(N + 2)(N + A) k 2 



D ^ ( ^ + 2)fee " A(Z) S^ )+ 4(, ElN + l) 2 

n n 

- S(zi) m 2 -Y, S(zi) d% - e A ^ fc( £ S^gn^d* 



i=i 



+8NaKke A(z) 



i=i 



TV 



8iV(iV + 2)aKA; 2 <5(«i)5(zj)^(«) =0, 



(20) 



where £ = 2N 2 {N + l)ank 2 . The modified Einstein equations in the bulk amount to giving a solution 
for k 2 : 



1-£ = =r/1 + 



8N(N + 1) a k 2 A 
(W + 2)(iV + 3) 



(21) 



An interesting special case arises for k = / \J2 N 2 ( N + Tjan (and hence £ = 1), where the two 
branches of the bulk solution (21) coincide. Of course, one require a > to ensure k > 0. Moreover, 
a fine tuning condition between the bulk cosmological term and the GB coupling a, 

(iV + 2)(iV + 3) 



8N(N + 1) a k 2 
7 



(22) 



is useful to study the localization of gravity to a solitonic 3-brane [7]. For N = 2, the choice £ = 1 
allows one to study a solitonic 3-brane solution with the supersymmetric interpretation [7]. The 
readers will note that (1 — £) = implies trivial corrections to the graviton propagator in (4 + N) 
dimensions, since all sub-leading-order corrections involve this factor. But one still recovers the RS- 
type braneworld solutions in one lower dimensions (i.e., in 3 + N dimensions). In fact, the choice £ = 1 
helps us to study gravity confined to a solitonic brane. It is interesting that the solutions with £ = 1 
have structures which are similar to the original Randall-Sundrum braneworld solutions. 

4 Finiteness of Newton's constant 

A finite four-dimensional gravitational coupling (or the Planck mass) ensures that the ordinary four- 
dimensional graviton is present as a massless and normalizable zero-mode solution of the theory. In 
the RS three-brane model of co-dimension one, the four dimensional Newton constant is determined 
to be 

~ M 2 Pl ~ Ml r dz e~ 3A ^ [l + e - 8a K5 e 2A ^ U' 2 - A"\ ] (23) 
G N Jo 

where £ = Aan 5 e 2A A' 2 . The last term [e"^ 2 )^'] ^° has to be finite for G$ to be finite. This is finite 
for a solution of the type A(z) = log(|z|/L + 1), since A'(0 + ) = 1/L and the background is invariant 
under z — > — z symmetry. 

Next consider the general case where N is arbitrary. In this case, one can find the four-dimensional 
Planck mass by reading off the coefficients of — m 2 in the linearized expression, equation (19) or 
equation (20), and then integrating over the extra dimensions. The result is 

N 



M 2 Pl 



/+oo 
d N z e ^ N+2 ^ (1 - + 8Na K k e A ^ £ 5{z t ) 



-oo 

M 2+N 2" (1 + k -W 



(N + 1)! ■' • (24) 
We will shortly justify this result by integrating the Ad part of the metric, as in the original RS model. 
For this purpose, one may replace rj^ in equation (7) by g^}{x), insert the metric into the action (2), 
and finally integrate over the Z{ coordinates. The net result is that 

S eff = M Pl J d 4 x yj-gW R + a(R 2 - AR 2 V + i?^ Ap ) + fine-tuned terms , (25) 

where the effective 4d Planck mass Mpi and a, the Ad analog of the L>-dimensional GB coupling a, 
are defined by 



M 2 Pl 



/+oo . 
d N z e -(JV+2)AW 1 _ 2N(N + 1)(JV + 2)an k 2 
-oo 

N 

+8{N + l)ctKke A W J2 S (*i) 
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= M * +JV (]vTI)!( 1 +^ fc " 7V ' (26) 

and, 

a = a^— / d N ze~ NA W. (27) 

M 2 pl J-OO 

Observe that to ensure a positivity of the four-dimensional graviton coupling (i.e. G4 > 0) and to 
preserve unitarity (i.e. a definite positive-norm state) condition a with GB term, one must take 
(1 + £) > and a > 0. 

The four-dimensional massless graviton corresponding to a bound state with the wavefunction, up 
to a normalization factor, is 

YWd~e-^+ 2 )^)/ 2 . (28) 

For the massless graviton mode m 2 = 0, the second square bracket in equation (20) then gives a 
non-trivial term 

8N(N + 2) auk 2 ^5{zi)5{ Zj )^{z) . (29) 

This exactly cancels with the last term in equation (20). Thus, the massless mode of the graviton 
fluctuation is unaffected by a Gauss-Bonnet action. And for \z\ » 1/k one always recovers the 
iV-dimensional Randall-Sundrum-type volcano potential. This, in turn, implies that the gravity in the 
bulk is not delocalized by a Gauss-Bonnet term even if N > 1. 

5 Graviton propagator in (4 + N) dimensions 

In order to study the graviton propagator, let us set N = 1 in (19) and momentarily change the bulk 
coordinate to z = e k ^/k. Then, in terms of A(y)(= k\y\), the Fourier modes e ip ' x ip(x,y) of the tensor 
fluctuation h^ u in five dimensions satisfy 

"(1 - 4a K5 A 2 ) (e 2 ^ I m 2 + d 2 - AA 2 + 2A) 

+4an 5 A e 2fc|s/| m 2 - 8aK 5 AA(d y + 2i)l ^(y) = . (30) 

Here dot represents differentiation with respect to y. To study graviton propagator for the matter- 
localized brane, one may place a test mass 27rm* on the brane z = 1/fc (or on the intersection of 
(N + 2)-branes) and ask for the corresponding Newtonian potential at a distant point on the brane. 
This may be done by inserting a source term G^ + n (27rm* )8 3 (x)5(z) on the rhs of equation (30). The 
five-dimensional graviton propagator Gs(x, z; x', z') is obtained by integrating over the Fourier modes 

g 5 (x,z;x',z') = J ^e^M^y). (31) 

We also need appropriate Neumann-type boundary conditions on the brane(s) to proceed further. 
For this purpose, we need to regularize the (5-function 1 . Also note that d y ip(0 + ) = —d y il>(0-), while 

1 This issue has also been explained in Ref . [20] , where a similar analysis is carried out by considering arbitrary order 
of curvature correction terms. 
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d y ip(\0\) = 0. Then the last term of equation (30)) gives 

- 1 -^^S(y)(d y + 2k)^(0 + ). 

Similarly, two other non-trivial terms arising from the first two brackets in equation (30), namely 
—AaKc,A 2 dyil)(y) and — 8aK 5 A 2 Aip(y), give 

-^aK 5 k 2 6(y) (d y + 2k)^(0+). 

The total non-trivial contribution at y = 0, from equation (30), is therefore 

2 (d y + 2k) V>(0+) - 8aK 5 k 2 (d y + 2k) ^(0+) + 8ok 5 k m 2 Y>(0+) . (32) 

This gives the correct Neumann-type boundary condition on the brane after integrating from just 
below to the just above the brane (i.e., in the neighborhood of y = 0). Hence a Neumann type 
boundary condition that one must impose on the brane at z = 1/k is 

e 

where £ = 4a' M~ 2 k 2 , and m 2 = q 2 = —p 2 . By the same token, one can derive (iV + 4) dimensional 
Neumann boundary condition at z = 1/k. As the calculation is much involved, we avoid here the 
technical details and simply give the final expression 



(zd z + 2 + ^-1- q 2 z 2 ) iP(z, z')\ z=1/k = , (33) 



^« + — +(]V>1) (T3£)J ^U/* = o. 



(34) 



This boundary condition was first introduced in Ref. [4], but for £ = case. One should also satisfy 
the following matching conditions at z = z': 

*P<\z=z' = i>>\z=z', 

d z - 4><) \z=z> = (i-^r 1 ^- (35) 

By satisfying (35) and (34) we can find the Neumann propagator in (4 + N) dimensions. To this 
end, we follow the derivations given in the Ref. [4, 14]. The formula for the Green function in five 
dimensions (N = 1), for both the arguments of the propagator on the brane, was given, for example, 
in [14], so here we give the expression valid for arbitrary N. The general expression for the Green 
function in (4 + N) dimensions for both the arguments on the brane is 



2(jv-i)(i+o 2 ( j (i+o 2 c! 



l-ZN + l 

(36) 



i + e w 

2 „„N 
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where q/k = w, c\ is a dimensional constant (c\ = 1,4, 64, for N = 1,3, 5, •••), but the last 
expression involving a logarithmic term will be absent for even numbers of extra dimensions (N = 
2,4, 6, • • •). We should mention that the second expression above will appear only for N > 1, though 
we have written it for arbitrary N. Note that only the first two terms do not vanish for £ = 1. 
We may split equation (36) into zero mode and KK mode 

,1,1. (N + l)k f d 3+N p e fr(*-*') 

-js.-H^-^- y (27r)3+jV g2 + (37) 



iv + (iv-2)e i 1 1 1 — e (Q/k) N ln / g . 



where 

(38) 

Equation (37) reflects a useful relation 

Gat+3 = - - — ~r Gn+4 (39) 

One observes from equation (37) that in the low-energy scale q/k(= w) « 1 and for N > 1 the 
contribution to the Newton potential from the continuum KK modes is almost negligible compared 
to the contribution of zero mode. One finds N = 1 as the simplest example. The limiting behavior of 
the Newtonian potential due to a point source of mass 2ir m* on the brane is therefore [14] 



27rm* f , „ ( 1 1\ 

, i 



'(5) 

G4 m* 



(40) 



(1 +£) 2k 2 r 2 

For £ = 0, this expression agrees with the result found in [3]. A difference of factor 1/2 to the leading- 
order potential correction from that of the original RS work [1] (which is the case of £ = 0) arises from 
the convention in a relation between four- and five- dimensional Newton constants. If one defines [1] 

K) = nr M ( 5 ) ( 41 ) 

instead of relation (39) the factor of 1/2 would be absent. To reconcile the result with an exact graviton 
propagator analysis, however, one may have to relate Gat+3 and Gjv+4 according to relation (39). For 
£ = 1, there is no correction from the continuum KK modes. Moreover, to avoid anti-gravity effect 
(G4 < 0), one requires £ > — 1, and another limit £ < 1 arises from the propagator analysis; this is 
consistent with the previous observation. Thus, for £ = 1 solution, there are no bulk propagation of 
graviton in AdS^+N space. By taking £ = 0, from (41), one recovers the relation = A; -1 and 
also recovers the Newtonian potential in a form identical to that obtained by Randall and Sundrum [1]. 
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For r » 1/k, the continuum KK modes with m » k, which are normally suppressed on the brane, 
would only sub-dominantly contribute to the potential generated by the 4d graviton bound state. 
For distances r « 1/k or at the origin, the continuum modes with m » k have unsuppressed 
wavefunction, thus as anticipated we get the (4 + iV)-dimensional potential even if £ > 0. 

A remark is in order. For £ = 1 and N > 1, only the first term in (38) can survive, but this is neg- 
ligible as compared to the zero-mode contribution. One, therefore, always recovers a localized gravity 
on the brane which precisely follows the Newton's law, and there may almost be trivial correction to 
the gravitational potential. We will show in next section that £ = 1 is indeed a viable solution of the 
RS-type braneworld models. 

6 Effective solutions with £ = 1 

It is not difficult to show that £ = 1 explains the RS warped braneworld compactification. It looks 
reasonable first to solve the modified Einstein equations in five dimensions. To this end, we shall 
introduce a brane action corresponding to a RS singular 3-brane with a positive brane tension A > 0. 
Thus the low-energy effective action in five dimensions can be taken to be 



S = Jd 5 x (J± - A + a (R 2 - 4Rp q R pg + R pqrs R™ rs )^ + j d A x y^J (-A) . 



(42) 



The RS single brane action is recovered for a = 0. The modified Einstein field equations following 
from (42) would simplify to give 

(A' 2 + A") (1-0 = ^k 5 5(z) (43) 

A' 2 -^A' 2 = -^Ae~ 2A ^ (44) 
2 12 y J 

where £ = 2ee 2A ^A' 2 , and e = 2aK$. One might have already noted that appealing to £ = in the 

bulk essentially implies A" + A' 2 = 0. In fact, £ = kills another possible solution, i.e. £ = 1, which 

indeed gives a physical solution. As is known, any braneworld compactification with £ = suffers 

from the so-called c-theorem [22, 23] based on Einstein theory. We should mention that the c-theorem 

may not be available to £ = 1 solution. We will explain this below only briefly, and detail discussions 

on this issue will appear elsewhere. 

One clearly sees from (43) that £ = 2eA' 2 e 2A ^ = 1 is a bulk solution of the field equations. From 

this we find 

e A( Z ) = r i dz 

Jo \[2e 



^ A(z) = rr^^M^' ( 45 ) 

where we have normalized the solution so that A(0) = 0. Since we have a non-trivial brane tension 
A > at the brane z = 0, equation (43) implies that we must satisfy the following relation on the 
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brane z = 0: 



\25{z) ( 2e 2 



l-H 8 ^) ' = Ak 5<5(z). (46) 



This, after regularizing the (^-function, uniquely determines the 3-brane tension, 



x = 2ikrr < 47 > 



Next, equation (44) can be used to fix the bulk cosmological constant 

3 



(48) 



8ttG 5 I 2 

The readers can easily check that £ = 1 consistently gives a physical solution even if the number of 
extra dimensions N > 1. It is interesting that e~ A ^ will converge as z — > oo, this is what we need 
for the solution to be physical. If one interprets the singularity at A" + A' 2 as the infrared cut-off (or 
the RS singular brane), then we may interpret the z — > oo, £ = 1 as the ultraviolet cutoff or anti-de 
Sitter boundary where gravity is strong. Changing to y-coordinate dy = e~ A ^ dz, the weak energy 
condition reads A"{y) > 0, and the warp factor reads A{y) = ±\y\/l. Since A'(y) takes a positive 
(negative) value for y — > +cx> (— oo), the no-go theorem [22] may not apply to the a > case. 

7 Completely localised gravity with £ = 1 

As we already emphasized by considering £ = 1, we can obtain equivalent N copies of the RS-type 
non-relativistic Schrodinger equations, but in one lower spacetime dimensions. We shall express the 
metric fluctuations again in a canonical form by changing the variables in equation (17) as h^ u = 
e (N+i)A(zi)/2 ^^ z ^ e ip-x e ^ Then the non-trivial part of linearized equations takes the following form 

N r / N 



I\ r / iV SI 

i=l L V i^j i^j 

(N- 1)(N+ l)(N + 3)k 2 \ 7 . . f%T _ . / (AT + 1) , \ 7 

+ 4^-^ttiF ) + {N + 1) agn ™ ( s « + S-^ ^ 



(49) 



Equation (49) clearly implies that the so called 'solitonic' solutions defined for £ = 1 have the structure 
which is formally the same as the "normal" or RS-type brane solutions, other than the last term. 
Some of the clear differences here are that (i) there is no bulk propagation of the graviton in (4 + N)- 
dimension, (ii) the solitonic brane is 5-function like, which is also seen from the overall coefficient 
in (49). The last term in equation (49) implies an appropriate jump condition(s) to be satisfied across 
the branes at Zi = 0, but does not contribute for z > 0. Thus, for £ = 1, the (N — l)-dimensional 
Schrodinger equation can also be written in the form 

4$6(zi) 



N(N + l)k 



^ 4 2 



ip = m 2 ip , (50) 
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where A = log (k ^ p \z p \ + 1) and i, j = 1,2 ■ ■ ■ , N, i ^ j. Thus, the bulk equations along Zj -directions 
are 



-dl + 



(N- 1)(N+ l)(iV + 3) 



■0 = m ip . 



(51) 



iV(iV + l)fc [ ~ Zj ' 4(\ Zj \ + l/k) 2 
The simplest choice is N = 2. Then (51) gives two five dimensional bulk equations of graviton, for 
each i,j = 1,2 but i / j. For the continuum of KK modes propagating along the solitonic 4-branes 
located at the Z\ and Z2 axes, the continuum mode solutions are given by a linear combination 



~ Vd^-I + V*0 [^Mi^-i + i/fc)) + i$)y 2 (m(|* j | + 

L,2. The coefficients 1?™ and Nm are determined fron 
(5-function normalization condition. These are readily evaluated to be 



(52) 

(i) (i) 

where j = 1,2. The coefficients Bm and Nm are determined from the boundary conditions and 



B U) = y i( m / fc ) N U) 

m ' J^m/k) ' m 



m 



4(i + ^ )2 ; 



(53) 



We consider the contribution of KK graviton exchange to the non-relativistic gravitational potential 
between two point masses mi, m.2 placed on the intersection of two orthogonal solitonic 4-branes with 
a separation of r. A total contribution of the Newtonian potential is therefore 



V(r) 
m\ rri2 

V{r) 



r 2 



G 4 



mi nT-2 



1 + 



k 2 r 2 



(54) 



where G4 = (3/4) x k 2 M~ A is used. In arriving at the second line it was important to use |^ m (0)| 2 ~ 
m j (4k) , which one can easily find from plane wave normalization conditions of the Bessel functions 
as in [5]. 

For N = 3, equation (51) would rise to give the following three six-dimensional bulk equations of 
the graviton, each bulk mode propagating along the solitonic 5-branes located at u(= z v — z q ) and 
v(= z p + z q ) axes. Hence 



2S( Zi ) 
3k 



-d 2 -d 2 + - 

u v (\v\ + l/k) 2 



ip(u,v) = - m 2, 4)(u,v) 



(55) 



where i = 1,2, 3, and i / p / q. In order to satisfy the appropriate boundary condition(s) implied by 
the 5-function potential at z,- L = 0, we must choose the linear combinations 



<Pm(u) 

Vm{v) 



Ni (m) \fu [ sm(m u u) + A m cos(m u u)] 



N 2 (m) y/(\v\ + l/k) J b/2 {m v {\v\ + l/k)) + B m Y 5/2 (m v (\v\ + l/k)) 



(56) 



(57) 



These solutions must satisfy the Neumann type boundary conditions at the brane junction u = 0, 
v = 0: 



u ■ 



d v -d u , 5 



+ 77 V' = 0, 



d v + d u 5\ 



(58) 
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For m 2 , ml > 0, satisfying the boundary conditions, the coefficients A m and B m are determined to 
be 



A mu = cot(m u u) , B mv = - 



The continuum wavefunction therefore reads 



VVn(j) = <Pm(Vj) X (p m (Vj). 



(59) 



(60) 



By considering the contribution of KK graviton exchange between two unit point masses mi, rri2, 
placed on the intersection of three intersecting 5-branes with a separation of r, one may find the total 
contribution of the Newtonian potential. This is given by 

V(r) n '■ QJ - 2 



mi ni2 
V(r) 



— + -rr M~ 5 / dm u dm v e - |^)(0)| s 



-G 4 



mi m2 



1 + 



(kr)- : 



(61) 



where yf m\ + m^ = in, G4 = (3/2) x k 3 M~ 5 and C\ is some constant of order 1. Again in arriving at 
this result, we have used the value |^ m (0)| 2 ~ m 2 /(6/c 2 ) determined from the plane wave normalization 
conditions of the Bessel functions, and specialized to the case where m u would extend down to m u = 0, 
and ^ mu (0) = 1. 



8 Comments on scalar and vector modes 



So far we have analysed the tensor mode h^ u of the metric perturbations, leading to graviton propa- 
gators and corrections to Newton potential. This analysis was important because a high-dimensional 
graviton may represent a continuum of four-dimensional states and the gravitational potential on the 
brane is mediated by an effective four-dimensional graviton. It may be relevant to know whether 
the scalar/vector modes of the metric perturbations have any significant roles in the higher dimen- 
sional brane background. One of the motivation to look after these modes is that in the conventional 
Kaluza-Klein theory of compact extra spaces, because of the isometries of a factorizable geometry, 
scalar (vector) modes correspond to the physical gravi-scalar (gravi-photon) modes [28]. But this does 
not seem to be the case for a non-factorizable geometry. 

The scalar and vector modes of the metric fluctuations, without any gauge choice, satisfy 



(1 - [d^d* (V " Vvwh) + (N + 2) d z A (d z h - 2d x h Xz + (N + 3) d z A h zz ) 

(1 - [d z (a A v - d„h) — (n + 2)d z Ad^h zz + d x d^h x]z 



= 0. 
= 0. 



Here h = h\, d z A 



k Eili s 9n{zi), and d z = d z 



j d ZN . 



(62) 
(63) 

(These 



1 and in the axial gauge, i.e., h z 







flZ j 



A'd z \z\ = k{20(z) - 1) 
equations were reported in [12, 14] but for TV 
Equations (62, 63) are trivially satisfied for £ = 1 solution. There are no terms having second 
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derivatives of the warp factor A(z), and hence there are no delta-function sources for scalar/vector 
modes. This is not the case for the tensor mode (see, for example, equation (20)), even after the 
setting £ = 1. Thus the scalar and vector modes of the metric fluctuations alone do not appear as 
dynamical fields of the theory [1, 7, 12], and hence may not be localized on the brane. This result is 
consistent with the observations made in [7, 21]. 

Since a massless graviton on the brane is assumed to be transverse and tracefree, it is reasonable 
to impose the gauge = = V/i^ to simplify the expressions (62), (63). For N = 1, since Z2 
symmetry on the brane already fixes h\ z = 0, one finds h zz = 0, thus the graviscalar and graviphoton 
components are pure gauge degrees of freedom. If the brane gravity is coupled to brane matter, one 
may still gauge them away using the brane-bending formalism discussed in [3, 4, 14] 2 . 

One can, of course, couple the scalar modes of the metric fluctuations to the scalar fluctuation 
of a non-trivial bulk scalar field. But the effective theory will be of scalar-tensor nature rather 
than conventional tensor gravity. Unlike the tensor modes, the scalar/vector modes of the metric 
perturbations are not automatically invariant in a gauge invariant manner. Indeed, for iV > 1, one 
may not set all the scalar /vector modes to zero or gauge them away, and in principle one can analyse 
them using some gauge-invariant formalism. The analysis is relatively simpler in a RS-type warped 
brane background with one extra dimension (i.e. N = 1) [21, 31], but the diagonalisation of linearized 
fluctuations of all metric fields for N > 1 is too hard, if not impossible. Some progress can be made 
by taking the simplest choice, N = 2, and using the gauge invariant formalism developed in [30]. We 
take the viewpoint that consideration of scalar and vector components would not change the results, 
at least the momentum structures of the graviton propagators, of this paper. 

9 Conclusion 

We have studied gravity localization to the RS branes for a more general class of N intersecting (2+iV) 

brane configurations by including higher order curvature terms in a Gauss-Bonnet combination. For 

N = 2, a GB term can support a feasible singularity at the brane junction, thereby allowing a 

non-trivial brane tension at the four-dimensional brane intersection. But this would not have been 

possible if the braneworld action had a contribution only from the Einstein term and the cosmological 

constant term [2]. The effective four-dimensional Planck scale Mpi has been derived, and also the 

explicit expressions for the graviton propagator in (4 + N) spacetime dimensions are obtained. We 

have explained in some detail interesting features of localized gravity with a Gauss-Bonnet term for 

'normal' or RS-type branes supported by (5-function-like sources, and also so-called solitonic branes 

defined for the class of solution £ = 1. 

2 See, for example, [29] for related comments about the cancellation of the scalar mode by the 'brane-bending' mode 
in the framework of a five dimensional 3-brane model. These comments may not be directly applicable to the 3-brane 
defined as the common four dimensional intersection of higher dimensional branes. 
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The present analysis might reveal some important features of the localized RS type braneworld 
solutions or solitonic braneworld gravity created by non-dynamical sources (i.e. in the absence of 
fields other than gravity). We have shown that the RS single brane model with a Gauss-Bonnet 
term in the bulk correctly gives a massless graviton on the brane as for the RS model. We provided 
a complementary description of the Newton potential corrections to the long-range interaction due 
to continuum modes living in the AdS space. Basically, our work extends the previous work in the 
literature [2, 4, 14] to the general (4+N) dimensional anti-de Sitter spacetimes, but the results reported 
here also include the contribution of the Gauss-Bonnet interaction term. The latter is the only viable 
combination of higher curvature corrections that one can introduce in the RS braneworld scenario of 
warped extra dimensions. 

For the class of solution £ = 1, the scalar and vector fields do not appear as the propagating degrees 
of freedom in the transverse z direction, though the zero mode of the tensor fluctuation is localised on 
the brane. This reveals that the physical brane could be a solitonic 3-brane living in six-dimensional 
spacetimes, rather than in five-dimensional AdS space. The number of non-compact extra dimensions, 
N = 2, is itself more interesting in the scenario of [2, 10]. 

In nutshell, the physical relevance of our analysis is two-fold - it is a meaningful extension of the 
RS model with an extra non-compact dimension to an arbitrary N, and - a natural generalization of 
the ADDK [2] scenario of infinite large new dimensions including the contribution of a Gauss-Bonnet 
term into the effective braneworld (gravitational) action. The present analysis also provides hints for 
a possible creation of solitonic 3-brane in a six dimensional warped bulk background, which would not 
be possible without the higher curvature terms. It is plausible that the solitonic braneworld solution 
defined with £ = 1 in N > 1 has a structure which is formally the same as for 'normal' or RS-type 
solutions. 
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